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Abstract 

The correct treatment of the close-to-boson character of excitons is known to be 
a major problem. In a previous work, we have proposed a "commutation technique" 
to include this close-to-boson character in their interactions. We here extend this 
technique to excitons with spin degrees of freedom as they are of crucial importance 
for many physical effects. Although the exciton total angular momentum may ap- 
pear rather appealing at first, we show that the electron and hole angular momenta 
are much more appropriate when dealing with scattering processes. As an applica- 
tion of this commutation technique to a specific problem, we reconsider a previous 
calculation of the exciton-exciton scattering rate and show that the proposed quan- 
tity is intrinsically incorrect for fundamental reasons linked to the fermionic nature 
of the excitons. 
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Non-linear effects in the optical properties of semiconductors have received consider- 
able interest, both experimentally and theoretically. In these non-linearities, the interac- 
tions between carriers play a crucial role. Up to now, two types of theoretical methods 
have been proposed to deal with these interactions. 

One method starts with the exact semiconductor hamiltonian written in terms of the 
free electron and free hole fermion operators and b^, and ends with the so-called semi- 
conductor Bloch equations or better, with an elaborate hierarchy of 2n-point density 
matrices ( 2 ~ 5 ^ which are then dynamically truncated to a given order in the radiation 
field. These approaches basically lead to a set of coupled equations for the time evolution 
of some expectation values of these free electron and free hole fermion operators. As the 
Bloch equations originally rely on the Hartree Fock approximation, various extensions 
have been proposed to include some correlation effects. In particular, it is possible to 
recover the low excitation regime, in which the free electron-hole pairs give rise to exciton 
resonances, by dropping all non-linear terms denoted as exchange processes ^\ However, 
being written in terms of free electron-hole pairs, these procedures are obviously appro- 
priate to systems well represented by free pairs. Even if various correlation effects can be 
included, these methods are mainly suitable at high density, when the screening is such 
that the exact correlations making the excitons are not crucial. For dilute electron-hole 
plasma however, these approaches can appear as starting from the "wrong" side, in the 
sense that, the excitons being the proper states at low density, they not immediately 
appear as the zero order terms. It is a priori rather subtle to extend these procedures 
with complete security to the low density regime in which the excitons play the dominant 
role. 

The other method seems much more appealing at first in the low density regime, since 
it relies on the fact that, in this limit, the electrons and the holes are bound into excitons ; 
as these excitons are made of two fermions, they very much look like bosons due to the 
spin statistics theorem. This idea leads to the so-called bosonic method in which the 
exact excitons are replaced by boson-excitons and the exact hamiltonian by an effective 
hamiltonian H c s in which appear boson-exciton operators. In addition to an obvious free 
exciton energy part, this effective hamiltonian must contain the interactions between ex- 
citons which are dominant at low density. Besides a "direct" term, which corresponds to 
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Coulomb interactions between excitons made with the same couples of electrons and holes, 
(e, h) and (e', h'), it has been shown ^ that the interaction part of H e s also contains an 
"exchange" Coulomb term which corresponds to Coulomb interactions between excitons 
made with different couples of electrons and holes, (e, h) and (e', h!) on one side, while 
(e, h!) and (e', h) on the other side. It is said that this exchange term has to be introduced 
in order to take into account the fermionic character of the excitons, i. e. the fact that the 
excitons are made of two fermions which can be coupled in different ways. We have been 
amazed to realize ® that this exchange term, quoted by everyone for 25 years, is incorrect 
since it induces a non-hermitian part in the effective exciton hamiltonian. Even if this dra- 
matic error is fixed, which is always possible at least intuitively, the major trouble with this 
bosonic approach actually comes from the fact that it misses purely fermionic terms ^ : 
Indeed the fermionic character of the exciton does not appear through Coulomb exchange 
terms only. 

This fermionic aspect is indeed very subtle, and it is necessary to use a "full-proof" 
procedure to derive all its consequences properly. The previous approaches, which all ^ 10 ^ 
end with incorrect results, were too intuitive to pick up all the consequences of this 
fermionic aspect properly. 

In short, we can say that : 

(i) On one hand, the semiconductor Bloch equations and their extention through 2n- 
point density matrices use "clean" free electron and free hole fermion operators. They 
have however problems with Coulomb interaction and its resulting correlations, which 
suffer from truncation procedures (even if they are somehow "controlled"). Indeed this 
Coulomb interaction must be included exactly in order to properly handle the evolution 
of the exciton bound states, i. e. the poles in the response function, which dominate at 
small density. 

(ii) On the other hand, the bosonic method introduces the exact exciton states a priori, 
but uses "too clean" boson operators to represent these excitons, so that it misses some 
important consequences of their underlying fermionic structure. 

What we would really like to do is to work with exact excitons a priori, i. e. electron- 
hole pairs in which the Coulomb interaction responsible for the bound states is put exactly, 
while we properly handle the fact that these excitons are not real bosons. This is this 
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approach we are now proposing. It has already been briefly reported in ref. (8), without 
the spin degrees of freedom for simplicity. 
The paper is organized as follows : 

In part I, we recall some very basic results on excitons in semiconductors to settle the 
notations. 

In part II, we introduce the Coulomb creation operator V} between the exciton % and 
the rest of the system and show how we can get out of it a direct Coulomb scattering 

tdir 
>>lnij ' 

In part III, we introduce the boson-departure operator Dy and the bare exchange 
coefficient \i ni j associated to it. 

In part IV, we show how we can produce various Coulomb exchange scatterings out 

Of Ztnij and X lnij- 

In part V, we use our commutation technique to calculate the matrix elements of the 
exact semiconductor hamiltonian between two-exact-exciton states. 

In part VI, we reconsider a previous approach ( n ) to the exciton-exciton scattering 
rate in which enter these H matrix elements and show why it cannot be correct. 

1 Exact semiconductor hamiltonian and one-exciton 
states 

The semiconductor conduction electrons correspond to S states with a s = ±1/2 spin 
along an arbitrary z direction. Let us call s the creation operator for such a free 
conduction electron with momentum k and spin s. For the free valence electrons, the 
situation is much more complex : According to the Kohn-Luttinger representation ( - 12 \ 
the valence energy matrix ( 13 ) is given by 

= Ak 2 + £?(k.I 3 / 2 ) 2 , (1) 

if one neglects warping, I 3 / 2 being the 3/2 angular momentum matrix. The diagonal- 
ization, which is done by taking the z axis along the k direction, generates the so-called 
light and heavy valence electrons. If we now add Coulomb interaction, we have shown ( 14 ) 
that this Coulomb interaction is not diagonal between heavy and light valence electrons, 
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a fact which seriously complicates the problem since it mixes the m = ±3/2 with the 
m = ±1/2 valence states. In quantum wells, one can forget about this complexity be- 
cause the ±3/2 and ±1/2 valence states having different hole masses, they are shifted 
differently by the confinement, so that the heavy-light hole Coulomb coupling usually 
gives a negligible contribution due to the splitting energy induced by the confinement, 
which appears in denominators. In bulk material however, the only consistent way to 
forget about these heavy-light hole Coulomb couplings is to assume B = 0, which corre- 
sponds to take all the valence electrons with the same mass. In doing so, is diagonal 
for any k direction, so that the valence electron states are simply characterized by a quan- 
tum number m = ±3/2, ±1/2, along a direction which is now arbitrary. Let us call b\^ m 
the creation operator for such a free hole with momentum k and angular momentum m. 
The exact electron-hole semiconductor hamiltonian then reads 

H = H + V ch + V cc + V hh , (2) 

where 

#o = E(A±4)4, s a k , s ±£e£O k , m , 4 = |^ e £ = S-> ( 3 ) 

k,s k,m e " 

A being the band gap. For quantum wells, m = ±3/2, while for bulk materials 
m = ±3/2, ±1/2, as a unique hole mass m h is then assumed for the two hole bands. 
Within this approximation, the Coulomb scattering of a hole is diagonal with respect to 
its "spin" m, so that a hole (k, m) is scattered into (k ± q, m) with the same m. The 
hole-hole, electron-hole and electron-electron interactions thus read 

^hh = l ~ E ^q 6 k+q,m 6 k'-q,m' 6 k',m'&k,m , (4) 

q^O,k,k',TO,m' 

V e h = - a k+q,s K'-q,m &k',m dk,s , (5) 

q^O,k,k',s,TO 

V ee = - V <i a k+q,s a k'-q,s' fl k',s' Ok,s , (6) 

Z q^O,k,k', S ,s' 

where Vq = in 3D and K, = in 2D, with V being the sample volume and S the 
well area. 

From the free electron and free hole operators a ks and b^^, one can construct the 
exciton operators B\ as 

B i = B li,Qi,Si,mi = XK k *l^i) a Kf,s, & L h .m, ' ( 7 ) 



where we have set 

Tfl h 

= k t + a e Q t , Kf = -ki + a h Qi , a e>h = e f— . (8) 

m e + m h 

Qi and i/j are the i exciton center of mass momentum and relative motion quantum 
parameter : In 3D, z/j corresponds to the quantum numbers (nj,/j,TOj), while in 2D it 
corresponds to two quantum numbers only (n i? TOj). (k\x Ui ) is the i exciton relative motion 
wave function in k space. In a similar way, it is possible to write the free electron-hole 
pair operator in terms of exciton operators : 

4.,* b L, mi = WhK - a e k h ) Bt itMhtSitmi , (9) 

as checked by inserting Eq (9) into Eq (7). Let us stress that this sum contains the bound 
states as well as the diffusive states. 

The one-exciton state B>\\v), where \v) is the electron-hole vacuum state, is eigenstate 
of H. Using Eqs (2-9), and the following Schrodinger equation for the electron- hole 
relative motion in k space, 

(4 + e£) (k|^) - £ Vq (k + q\x v .) = e„.(k\x v J , (10) 

q^O 

where e Vi is the energy of the relative motion state \x Vi ), it is indeed easy to check that 

HB\\v) = (H + V ch )Bj\v) = Et B}\v), E t = A + e Vi + ^ . (11) 

2{m e + m h ) 

Let us again stress that in 3D this result strongly relies on the fact that the Coulomb 
interaction is diagonal between holes, i. e. it does not contain terms in &k/_ qim /frk',m with 
m! 7^ m. 

In many problems dealing with excitons, the photons play an important role. As a 
a + photon with spin (J = 1, M = ±1) creates an exciton with the same total angular 
momentum, it may appear appropriate to use, instead of the (sj,TOj) angular momentum 
variables of the separate electron and hole, the ( Jj, Mj) total angular momentum variables 
of the exciton, with ( Jj = 2, Mj = ±2, ±1, 0) or (Jj = 1, Mj = ±1,0). One can easily go 
from one set of operators to the other by 

B Lq M = T,(^^\^M l )Bl Qi ^ mt , (12) 

s i )Tn<i 

^.Qi.-i.mi = T,( J i> M i\ s i> m i) B ii,Qi,Ji,Mi- ( 13 ) 



Ji,Mi 
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Only two of these Bl q j M . states are coupled to light, namely (Jj = 1, Mj = ±1), 
while the six other exciton states correspond to the so-called "dark" excitons. However, 
even if these dark excitons are not coupled to light, they are generated by Coulomb 
scatterings so that we must keep these 8 operators B\ ^ j M . anyway in order to possibly 
describe the exciton-exciton scatterings properly. Moreover, as the Coulomb interaction 
is diagonal within the (s,m) quantum numbers, but not within the (J, M)'s, it is in fact 
far simpler to work with the 8 operators BIq for all processes dealing with Coulomb 
interactions, and just at the beginning and the end of the calculations, use Eqs (12-13) to 
transform these B^q operators into the bIq JM operators coupled to the light. This 
is why all over this work, we will use the B^q operators only. 



2 Coulomb creation operator V? and direct Coulomb 
scattering Qj^. 

In standard problems with interactions, one usually divides the system hamiltonian H 
into H + V, where H is the so-called "non-interacting" part, i. e. the part of H which 
can be diagonalized, while V is the interaction part, i. e. the part of H which cannot be 
handled exactly, but is hoped to be small enough to be treated as a perturbation. In 
addition, when H can be written as H + V, the H eigenstates form an orthogonal basis 
which can be used to expand any state of the system. 

For interacting excitons, even if we could guess that the energies Ei of individual exci- 
tons should appear in the "non-interacting" part of an hypothetical exciton hamiltonian, 
one cannot divide the electron-hole Coulomb interaction given in Eq (5), into a part which 
would bind a specific electron to a specific hole to form the % exciton, and a "rest" which 
would make this exciton to interact with other excitons : The electrons, as well as the 
holes, being indistinguishable particles, such a formal separation is indeed impossible. 

A separation, similar in its spirit to the splitting H = H + V, is nevertheless possible 
through our commutation technique. If H is the exact semiconductor hamiltonian, given 
in Eqs (2-6), and if B\ is the exact exciton creation operator, given in Eq (7), we find 
that their commutator reads 

= EiBt + V}, (14) 
7 



where is given by 

W -1 = E a l-<l,s a P,s ~ E 6 P-q,™ V™ • ( 15 ) 
p,s p,m 

The coefficient 7;i(q), characterizes the scattering of a z/j exciton into a ^ state under a 
q Coulomb excitation. Using the following relation, 

£>„|k + aq) (k|av) = (a^e^M , (16) 

k 

it is easy to show that it is given by 

lH {q) = {x Vl \e ia ^ T - e~ ia ^ r \x Vi ) . (17) 

The analogy with the usual separation if = H + V is quite transparent if we note 
that Eq (14) also reads 

HB\ = B\(H + E l ) + V? . (18) 
By considering the state Bj\<j)), where \<p) is any electron-hole state, Eq (18) leads to 

HB\\4>) = EiB\\<l>) + B\H\4>) + V t ] \4>) . 

In the first term, the contribution of the % exciton to the energy of the system is just Ei 
as if this exciton were not interacting with the other electrons or holes of The second 
term corresponds to H acting on \<f>) independently of the presence of the i exciton : The 
operator {H + on the right hand side of B\, thus plays the role of the H part of the 
hamiltonian for usual problems in which H can be written as H + V. The third term is 
there because the i exciton does in fact interact with \<p). This operator thus describes 
all possible Coulomb interactions between the % exciton and the rest of the system. It has 
to be seen as the formal equivalent of V in H = H + V. It is however important to stress 
that, while the usual V's conserve particles, i. e. contain the same number of creation 
operators a) and destruction operators a, the operator V} is not a real potential in the 
sense that it contains one additional electron-hole pair creation operator B^ (see Eq (15)). 
This is why we call it Coulomb creation operator. In addition, usual two-body potentials 
have prefactors which depend on four indices which are the ones of two initial and two 
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final states. We here see that the prefactors appearing in V} depend on two indices only, 
I and i. 

In order to cope with these difficulties and generate a "scattering" which depends 
on four indices, we can push the commutation technique one step further and calculate 



Using Eqs (7-9) and (15), we find 

V i> B j]= E ^q7K(q)7nj(-q)^ )Q . +q)8 . )f ^sI n)Q ._ q>s . )ro . , (19) 

q^O,i/j,i/„ 

which formally reads 

V i> B j] =Y,ZtniJ B l B l ■ (20) 

l,n 

If we compare Eq (19) to Eq (20), and symmetrize the result with respect to (l,n) (which 
will appear convenient afterwards), we find that the direct Coulomb scattering can 
be written as 

& = & = & = \ [A dir fe) ^ fa) + (/ - n)} . (21) 

It contains an angular momentum contribution and an orbital contribution which appear 
as two independent factors. For direct processes, this angular momentum contribution, 

A dir fl i \ _ r r r r f22) 

\ n J/ s i> s i u mi,mi u s n ,Sj u m n ,mj i K^"^) 

just says that the I exciton has the same electron and hole momenta as the i exciton 
and similarly for the n and j excitons. Let us stress that this angular momentum part 
would be much more complicated if the (J, M) variables for the exciton total angular 
momentum were used ; this is why it is indeed appropriate to keep these (s, m) variables 
as long as we deal with Coulomb scatterings even if they are not the good variables for 
semiconductor-photon interaction. 

The orbital part of this direct Coulomb scattering is given by 

i dlT (!nj) = S Ql+ Q n ,Q i+Qj VQ ( _ Qi 7K (Q, - QO 7 n,(Qn " Q,) • (23) 

It of course contains the fact that the center of mass momenta have to be conserved, 
Qi + Qj — Qi + Qn, m the scattering. It also contains the two factors which characterize 
the scatterings of one exciton from a state to a v\ state and the other exciton from a 
Uj state to a u n state, under the Coulomb excitation Q; — Qj = — (Q n — Q^). 
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The appendix A contains the explicit calculation of these 7^(q) factors for the lowest 
S states. 

Even if the above expression of £ dir (n*,-) * s Quite transparent and extremely convenient 
for explicit calculations, it will appear useful, when we will generate exchange Coulomb 
processes, to note that this direct Coulomb scattering is also equal to 



£ dir W = E W k *> MK) O^M <b\x VJ ) 



X E [^Kf.Kf+q S K h >K h - (e <-> /l)] [<JK«,KS-q ^Kft.K* ~ (e <"> /l)] • (24) 
q^O 

In appendix B, we show that this £ dir {^ n l ^j is nothing but the direct part of the 
Coulomb interaction widely quoted ( 7 ' 15 ) in the effective boson exciton hamiltonian, namely 



I dir (J» i j ) = J dr e dr h dr e , dr h , <j>*(r e , r h ) fafe, r h .) 

x [V ee , + V hh , - V eh , - V e , h ] 0i(r e , r h ) (f)j(r e/ , r h ,) 



cdir ji I \ 



(25) 



where 0j(r e , r h ) is the whole wave function of the i = (z/j, Qj) exciton, 

(26) 

we have found a formal 



(r e ,r h \B}\v) = &(r e ,r ft ) = _L e * Q -^+^)(r e - r,|^) . 



and 



Bt' 



By these two successive commutators H, B\ 
way to generate a scattering of two excitons into two other excitons (l,n). If we 
look at its expression given in Eq (25), we see that the coefficient £ dir (^) obtained 
by this procedure turns out to have a very simple physical meaning : It corresponds to 
all electron-electron, hole-hole and electron-hole Coulomb interactions between and 
(/, n), when the initial and final excitons are made with the same electron-hole pairs (e, h) 
and (e', h'). 

The scattering processes corresponding to £ dir fyj) are shown in Fig. (1) : One of the 
electrons, or holes, suffers a (Qi — Q«) Coulomb excitation, while another particle has a 
— (Qz — Qi) excitation, the excitons before and after the scattering being made with the 
same electron-hole pairs. This is why we call it a direct Coulomb scattering, by contrast 
with the exchange Coulomb scatterings which will appear below, in which the excitons 
before and after scattering are made with different pairs. 
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3 Boson departure operator Di 3 and bare exchange 
coefficient A/ n ^ 

Even if the above procedure allows to formally extract from the whole electron-hole 
interaction the part which corresponds to a Coulomb interaction between excitons, such 
a formal extraction seeming not obvious at first, we are far from having picked all the 
physics which controls the interactions between excitons. Another rather subtle origin of 
these interactions comes from the fermionic character of the excitons. From handwaving 
arguments, we can say that two excitons feel each other even in the absence of any 
Coulomb interaction, because they are made of two fermions and the fermions of two 
excitons must be in different states. This condition in itself produces a "link" between 
excitons. There is no need of Coulomb forces. It is thus very likely that the exciton- 
exciton interaction has to contain a purely fermionic contribution, quite different from 
the Coulomb interaction dressed by exchange processes as thought by everyone up to 
now. 

In this quite tricky determination of the correct exciton-exciton interaction, it is how- 
ever highly necessary to formalize the above handwaving argument. We want to find a 
coefficient, which depends on four exciton indices (Inij), and which originates from the 
fermionic character of the excitons only. This coefficient has of course to be linked to the 
fact that the excitons are not real bosons. The most direct way to have this property 
appearing, is to start with the commutator of two (exact) exciton creation operators. It 
leads to the boson departure operator Dij defined as 



Bi,B}] . (27) 



Using Eq (7) and the standard commutation rules for the fermion operators a' and b', we 
find that this operator reads 



kj ,k 



D ij = H (zi/Jki) (kj\x Vj ) S mi>mj <W K h a K e )S . a K -, Sl + S Si>Sj <?k?,k| b Kh b K h m . 

3 3 j i ' 3 



• (28) 



= if we have both Sj ^ Sj and m 8 ^ rrij. 

From this boson departure operator, we can construct a four index coefficient by taking 
the commutator of this D^ operator with an exciton operator. Using Eqs (7,9,28), we 
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find 

D lh B}} =2j2^im 3 Bl , (29) 

n 

where the bare exchange coefficient \i ni j can here again be split into an angular momentum 
part and an orbital part : 

1 r 



^Inij ^nlij ^Inji 



. cxch 



fr^) A fr^) + (I <- n)] . (30) 
The angular momentum part, 

^ fr j) = ^ s h s i ^ m h m j ^s n ,Sj 3m n ,mi , (31) 

just says that the spins of the / and % electrons are the same while the angular momentum 
of the I hole is not the one of the % hole, as in A dir fr j), but the one of the j hole. The 
orbital part appears as 

A (r/i) = Z M k l) (^nl k n> (ki^i/i) (kj|^) fef.K? £ Kj h K * 5 K-,K e <5 K jj )K h . (32) 

In appendix B, we show that this orbital part A fr^ can be rewritten in r space as 

A fr^) = J dr e dr h dr e , dry $"(r e , r h ) 0*(r e /, r h >) &(r e , r h >) <pj{r e >, r h ) = [A (} '„)] . (33) 

From the above expression, we see that, in this orbital part, the electrons and the holes 
forming the excitons are exchanged, the / and % excitons being made with the same electron 
but different holes. Note that the (I <-> n) term of Eq (30) restores the (e, h) symmetry 
in Xi n ij, which is broken in A fry . 

Eq (33) shows in a transparent way that this A fry exchange coefficient can exist 
because of the composite nature of the excitons which can be formed in different ways, 
the electrons and the holes being indistinguishable fermions. An equally transparent 
link between this coefficient and the composite nature of the excitons can be obtained 
by considering the two-exciton operator B\ Sj, each being given by Eq (7), and by 
binding the electrons and holes of these two B^s in a different way through Eq (9). This 
gives 

B\ B) = - £ \ lni3 Bj i£ , (34) 

l,n 

with the same prefactor A lni j as the one of Eq (29). 
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n j 



The electron-hole exchange appearing in the first term of Eq (30) is shown in Fig. (2). 
It corresponds to cross the holes when forming the or (l,n) excitons : It is indeed a 
bare exchange process in the sense that it does not contain any Coulomb interaction, by 
contrast with the Coulomb exchange processes which will appear later on. 

The four <5's appearing in A i^ n l j) given in Eq (32) impose the expected exciton total 
momentum conservation Qj + Qj = Qi + Q n as well as kj + kj = ki + k n . By expliciting 
these four <5's, we can rewrite A (li t j) m a compact form, 

X (j, ;) = S Ql+Q n ,Qi+Qj Flnij MQl ~ Qi), UhiQn ~ Qi)) , (35) 

where Fi ni j(p, p') is a sum over one momentum only, 
JWp,pO = E(^l k -^> KJk+^) (k+P^|^) (k-P^|^.) . (36) 

k 

Although less transparent for physical understanding, the above expression of A 
is more convenient for explicit calculations. When the four excitons are in S states, the 
function Fi ni j(p, p') a priori depends on three independent parameters, p, p' and the angle 
9 between p and p'. Some interesting values of these parameters are : 

(i) p = or p' = 0, which corresponds to Q/ = Qj or Q n = Qj : The final excitons 
have the same momenta as the initial excitons ; 

(ii) 9 = or 9 = ir, which corresponds to p parallel to p'. Q/ — Q« is parallel to 
Qn — Q« = Qj — Qi in the particular case of Q« = Qj or Qz = Q n , i. e. when the 
two initial or the two final excitons have the same momentum. In this case, the two 
independent parameters can be taken either as a e |Q; — Q?| and a^|Q n — Qi\, or better as 
the momentum transfer \Qi — Qj| and the mass ratio m e /mh- 

In appendix C, we have calculated Fi ni j(p, p') when all the z/'s are equal to the IS 
ground state and Qj = Qj = Qi = Q n - We get 

F lslslsls (0,0) = (33n/2)a 3 x /V in 3D 

= (An/5)a 2 JS in 2D , (37) 

a x being the 3D exciton Bohr radius. 

We have also numerically computed A (V-,) in the 2D case, when v\ = v n = z/j = Vj = 
IS, Qj = Q, = 0,Qi = -Qn = P- Fig. (3) shows 

MP) = T^A (4 °) , (38) 



2D 
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for three values of m e /mh, as a function of P = p\2D, with A2D = a x /2 being the 2D Bohr 
radius. We see that A(p) weakly depends on the mass ratio, and that it is significant for 
values of p such that pA 2 D < 4 only. 



4 Coulomb exchange scatterings 

From Eq (33), we see that the \i ni j coefficient is dimensionless. So that it has to be 
"cooked" with quantities homogeneous to an energy in order to possibly appear in an 
exciton-exciton scattering. From the two energy-like quantities we have yet found in this 
problem, namely the bare energies of the excitons and the direct Coulomb scattering 
fifnij, there are of course various ways to construct a \i n ij dependent scattering. In the 
next section, we will show that two specific combinations of the bare exchange coefficient 
\inij and the direct Coulomb scattering appear in a natural way, namely 



£lnij £lnpr ^prij j (39) 

p,r 

£lnij = E ^Inpr £ pr ij ■ (40) 
p,r 

They are shown in Fig. (4). 

The simplest way to calculate these sums is to use Eqs (30-32) for \i ni j and Eqs 
(21,22,24) for ■ We find 

= \ (A cxch (^) f isht (V,) + (I <- n)) , (41) 

where the angular momentum conservation part A cxch is just the one appearing in 
the bare exchange coefficient \i ni j (Eq (31)), while the orbital part of this right exchange 
Coulomb term is given by 

r sht (^) = Ee dir (^)M^) 

^ (avjk/) (x Un \\t n ) (\ti\x Ui ) (\tj\x Uj ) 



p,r 



X V <i ^Kf,K e +q<5 K f,Kj - ^Kf,K^+q^Kf,K?J [^K= ,KJ-q^K^K^ ~ ^Kh,Kj-q^K= ,K?J • (42) 

Note that the second term of each bracket does not correspond to (e <-> /i) as for £ dir (j/j) 
given in Eq (24). 
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The set of four <5's appearing in £ nght (^\) again impose the expected momentum 



n j 



conservation Qj + Qj = Q/ + Q„. By expliciting these four 5 functions, it is possible to 
rewrite £ right (^) in a more compact form in terms of sums over q and one k only. In 
addition to the function Fi ni j(p, p') appearing already in the bare exchange coefficient 
(Eq (36)), two other functions enter £ n s ht {^ n l ^j. They are defined by 

W ( nn '.„\-sru n> p + p\/„ iu , p + p\/u , P-P 7 , »\„ wi. P~P ; 



G^(p,p';q) = E(^Jk-^) (xjk+^) (k+^+q|x„> (k- = 



(43) 



^(P,p';q) = E(^Jk-P±^) (^Jk+P±^) (k + P-^K) (k-P-^ + ql^,) . 

(44) 

(We can see that Fj ni j(p, p') is also either G^(p,p';0) or G^ -(p, p'; 0)). In terms of 
these three functions, the right exchange Coulomb term reads 

|nght = 5 Qi+ Q n Q. + Q. ^ Vq 

q^O 

x [F Jnij (a e (Qz - Qj) - q, a h (Q n - Q*)) + F lnij (a e (Qi - Qj), a h (Q n - Q<) + q) 
-G% (a e (Q, - Qi),a h (Q n - QO; -q) - G^- (a e (Q, - QO, a h (Q n - QO; q)] -(45) 

This £ right coefficient can be rewritten in a form much more transparent for the 
physical understanding, although much less convenient for calculations. In appendix B, 
we do show that 

X \Veef + V hh > - V M - Vefh] 0j(r e , Y h ,) 0j(r e /, T h ) . (46) 

The physical meaning of this £ nght (f n l j) becomes now clear : It contains the electron- 
electron and hole-hole interactions between two excitons when these excitons are made 
with their electrons and holes coupled in a different way. With respect to the electron-hole 
interactions, the situation is however more subtle as these interactions are between the 
excitons on one side, but inside the excitons of the other side, due to the exchange of 
the electrons or holes making the excitons of the two sides. This £ right (^ n l j) turns out to 
be nothing but the exchange Coulomb interaction which appears in the exciton-exciton 
scattering coefficient of the effective boson exciton hamiltonian widely quoted ^ . 

The coefficient reads as Cinij* except that V e h' + V e 'h is replaced by V e h + V e 'h' in Eq 
(46), so that the electron-hole Coulomb interactions are now between the excitons 
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and inside the (/, n) excitons. As a consequence, the two coefficients £ left and £ n s ht verify 



i 



left 

Inij 



fright 



(47) 



with a similar relation between the orbital parts £ left and £ ri s ht . This can also be seen 
from Eqs (25,33,39,40). 

From Cinij 1 an d £]niji ^ i s possible to construct a quantity symmetrical with respect to 



the electron-hole interactions, by taking 

pcxch _ j_ 



fright 
SZnij 



tleft \ 



(48) 



Using Eq (46), the corresponding orbital part of this Coulomb exchange coefficient reads 



*cxch 



(ij) = J dr e dr h dr e > dr h > 0*(r e , v h ) 0*(r e /, r h >) 



x 



V ee < + V hh , - -(V eh + Vjh + V M + Vefh 



(r e ,r^)0j(r e ',r A ) . (49) 



It contains the four possible electron-hole Coulomb interactions between two electrons 
and two holes. This ^f^J 1 could be a reasonable scattering for the interaction part of an 
hypothetical effective exciton hamiltonian : As ^fj 1 = (^ijfn^) > its contribution would be 
hermitian, which is not the case neither for nor for alone. 



5 Matrix element of H between two-exciton states 

In some physical problems such as the one considered in section VI of this paper, we 
are led to consider the matrix elements of the exact hamiltonian H between two-exciton 
states, namely (v\BiB n H B\ Bj\v) . Using Eqs (14,20) and V^\v) = 0, as obvious from Eq 
(15), we find 

HB\B]\v) = {E i + E j )B\B)\v) + Y,C% B l B t\ v ) ■ (50) 

p,r 

So that this matrix element reads 

(vlB^HBjBjlv) = (Ei + E^B^BlB^v) + Y.tfriM B i B n B l B t\ v ) ■ ( 51 ) 

p,r 

To go further, we must calculate (v\BiB n B\ Bj\v) . This is easily done using Eqs (27) 
and (29). Since D^v) = 0, as obvious from Eq (28), we get 

(vlB^BjBjlv) = 5 H 5 nj + 5 ni 5 tj - 2X lnij . (52) 
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The two first terms are naive : They differ from zero for (l,n) = i. e. when the 

excitons on both sides are identical. The last term is more subtle. It has a fermionic 
origin : In addition to the fact that the coefficient \i ni j is directly related to the boson 
deviation operator D H through Eq (29), its fermionic origin can also be traced back to Eq 
(34) which says that a product of two exciton operators writes as a sum of products of any 
two other B^s so that B\B^ in fact contains a "piece" of any B\B\, even if (l,n) 7^ 
In other words, the B}Bj\v) states do not form an orthogonal basis for two-pair states. 
By inserting Eq (52) into Eq (51), we get 

(vlB^HBjBjlv) = (E, + Ej)(5 u S nj + S u 8 ni - 2\ lnij ) + 2(^, - £g) . (53) 

This calculation thus produces the Coulomb exchange term A similar calculation 

done with H acting on the left gives 

(vlB^HBjBjlv) = (E l + E n )(S li S nj + S lj S ni -2X lnij ) + 2(^ ij -^ t ) . (54) 

Here appears From these two expressions of the H matrix element, we deduce that 

and are linked by 



C -^ff = (E l + E n -E i - Ej)X lnij , (55) 

so that they are equal for Ei + E n = Ei + Ej only. 

It is possible to check that Eq (55) remains valid when £ and A are replaced by their 
orbital parts | and A. Using this equation, we can easily calculate the difference £ left — £ right 
and compare it to £ lcft (or £ n s ht ). Ciuti et al ( n ) and Rochat et al ( 9 ) have calculated 
£ left i^nj) for 2D scatterings between v { = Vj = IS, Qj = Qj = 0, and v x — v n — IS, 
Qi = — Q n = p. Fig. (4) of reference (9) precisely shows 



as a function of P = p\2D, with \ 2 d being the 2D Bohr radius. On our Fig. (5) is shown 
A9{P) = i^~ D P Cft °°) " ^ ^ °)1 = T W^Afc,) , (57) 



with A(p) given by Eq (38). We see that, except for very small pA 2 D, this difference is of 
the order of g cxch (p), unless m e « m h . 
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We can rewrite the H matrix element between two-exciton states in a more symmet- 
rical form with respect to (l,n) and by taking half the sum of Eqs (53) and (54). 
If in addition, we consider this H matrix element between two-exciton states normalized, 
but still not orthogonal, namely 

RtotLA 1 

\*l>n) = \A 3 ) = MM r = TT B ! B i\v) > (58) 

where Mij = (1 + <% — 2A ij -jj) 1 / 2 due to Eq (52), we find, using Eqs (53-54), 

^ ln \H\Aj) = 1 -(E l + E n + E l + Ej) 5 lnij + - Ct - (1 - 8 lnij )r ]lni3 ] . (59) 

5i ni j is equal to 1 when the excitons on both sides are identical i. e. (l,n) = and 

otherwise, while r]i ni j is given by 

Vinij = ^(E l + E n + E i + Ej)\ lni j . (60) 

The first term of this H matrix element just corresponds to the energies of the non- 
interacting or (l,n) excitons if, in addition, the two-exciton states are assumed 
to be orthogonal. The last terms of this H matrix element come from interactions. They 
have three origins. The first term, is a direct Coulomb scattering between excitons 
made on both sides with the same electrons and holes. The second term, £f*^ h , is an 
exchange Coulomb scattering in which the holes (or electrons) making the excitons are 
exchanged. The last term, r]i ni j, has a purely fermionic origin : It is directly linked to 
the fact that the \ipij) and \ipi n ) states are not orthogonal for (l,n) ^ It does not 

exist in diagonal scatterings. Let us however stress that, as i]i ni j depends on the sum of 
the four exciton energies, it is band gap dependent ; so that it is very unlikely that such 
a i]i ni j term appears in a physical scattering. We are going to come back to this problem 
in the next paragraph. 

It might be useful to mention that the diagonal matrix element of H reduces to 

<lM# llM = E l + E 3 + {^ J - C-')/(l + - 2X tJZJ ) , (61) 

with possibly replaced by ^f] 1 * or as these three scatterings are equal for diag- 
onal processes. (Ei + Ej) is just what we would expect for the expectation value of the 
hamiltonian between two non-interacting boson-excitons i and j. The other terms come 
from scattering processes between the % and j excitons resulting from Coulomb interaction 
and Pauli exclusion. 
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6 Exciton-exciton scattering rate 



Let us end this work by considering a physical quantity in which these H matrix ele- 
ments between two-exciton states may appear in a direct way. Ciuti et al (CSPQS) have 
proposed ( n ) to calculate the elastic Coulomb scattering of Is excitons by Is excitons, 
namely 

(Is, Q, S) + (Is, Q', S') - (Is, Q + q, S f ) + (Is, Q' - q, S' f ) , (62) 

(see their Eq (3)), from the matrix element of the exact hamiltonian H between these 
two-exciton states. Let us rewrite this Eq (62) as (i) + (J) — > (I) + (n). In the following 
we will drop the angular momentum parts for simplicity, as they are unimportant for the 
problem we raise. 

The one-exciton wave function 0Q(r e , r h ) given in their Eq (1) is just our 0«(r e , r^) (see 
Eq (26)), while the two-exciton wave function 0qq/ given in their Eq (4) reads, within 
our notations, 

(f)ij(r e , r h , rv, r h >) = - [&(r e , r h ) (f)j(r e ,, r h ,) + (r e , r h <-> r e ,, r v ) - (r e <-> r e ,) - (r h <-> r h /)] , 

(63) 

so that it is just the wave function of the state B}Bj\v). 

It is possible to show directly from Eq (63), i. e. their Eq (4), that 

{4>ln\4>ij} = Sli $nj + $lj $ni ~ ^Inij , (64) 

with 

knij = ^Jdr e dr h dr e , dr h , [0,*(r e , r h ) (f)* n (r e /, r h >) + (I <-»• n)] ^(r e , r fc #) fafe, r h ) , (65) 

which is nothing but our Eqs (30,33) and (52). As a consequence, the |0/ n ) states are 
not normalized. Eq (64) also shows that the initial and final states of the exciton-exciton 
scattering considered in Eq (62) are not orthogonal for q ^ or q ^ Q' — Q. So that, 
even if for such q, the S terms of Eq (64) give zero, this matrix element is not zero due 
to the \i n ij term. 

According to CSPQS, the scattering amplitude corresponding to the process of Eq 

s s f 

(62) should be equal to H S g, f (Q, Q', q) = {(j)ij\H\(j) ln ) (see their Eq (6)), where H is the 
exact two-electron and two-hole semiconductor hamiltonian. Their Eq (6) is nothing but 
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(v\BiBjHB}B^\v) written in r space. Our equations (53,54) immediately give 

(vlB^HBjBUv) = 2 [SX - - (E l + E n )\ jln ] (66) 
= -(Ei + Ej)^] , (67) 

as the 5 terms give zero for q ^ and q ^ Q' — Q. If we consider the results given 
by CSPQS in their Eq (7), we see that the first two terms given by their Eqs (12-14), 
exactly correspond to the two terms of (2£^[ n ) (see our Eqs (21,25)), while the last two 
terms (their Eq (16) and its "hole" equivalent) correspond to the two terms of (— 2^[jf^ t ) 
(see our Eqs (41,46)). Consequently they should have found the result given in Eq (67). 
Even if calculations in r space are quite cumbersome, it is in fact possible to check the 
existence of the missing term (Ei + Ej)\iji n directly from a (tedious) calculation of the 
right hand side of their Eq (6). The origin of their mistake probably comes from the fact 
that they have considered the two states \<f>i n ) and \4>ij) as orthogonal, which is not true 
(see Eq (64)). 

The correct value of the (4>ij\H\4>i n ) matrix element given in Eqs (66-67) leads us to 
question the validity of the whole procedure to determine the exciton-exciton scattering 
rate. 

(i) First the \(f>ij) wave functions are not normalized, so that it would be reasonable 
to use normalized states and relate the exciton-exciton scattering not to (0y|iJ|0/ n ) but 
to (0ij|-ff|0«n)/(0jj |0jj) 1//2 (0/n|0«n) 1//2 - This would add a prefactor to all scatterings (see 
Eqs (58-59)). 

(ii) There is however a much more dramatic problem with CSPQS procedure. The 
existence of the last terms of Eqs (66-67) is highly non physical for an exciton-exciton 
scattering : As it contains the exciton energy E t + E n (or E,- L + Ej), which is essentially 
equal to twice the band gap, such a scattering would be band gap dependent. We could 
of course get rid of this problem by deciding to use for H an hamiltonian without the 
band gap i. e. by dropping A in Eq (3). 

Actually the exciton-exciton scattering rate cannot be related to this matrix element 
for more fundamental reasons. In usual problems dealing with interactions, the hamilto- 
nian can be written as H = H + V. The Fermi golden rule then says that the transition 
rate between two different H eigenstates \i) and |/) results from the possible interaction 
between these two states through |(/|^|i)| 2 . When these states are eigenstates of H , 
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they form an orthogonal basis, so that this transition rate is also equal to \ (f\H\i)\ 2 , since 
(f\H \i) = for \f) 7^ \i). All the difficulties here come from the fact that H cannot be 
written as H + V : There is no hamiltonian for which the two-exciton states \4>i n ) and 
\4>ij) are the (exact) eigenstates, so that (i) these two-exciton states are not orthogonal 
for (l,n) 7^ (i,j) ; (ii) the matrix elements of H between these states have no reason to 
be equal to the matrix element of "the" interacting potential V - which, anyway, can- 
not be formally extracted from H. Consequently there is no reason to believe that the 
exciton-exciton scattering rate is given by \{<pi n \H\(f)ij)\ 2 . 

7 Conclusion 

We have reconsidered our commutation technique designed to deal with interacting close- 
to-boson particles and extended it to excitons with spin degrees of freedom. Although 
more cumbersome, the expressions of the two important parameters of this commutation 
technique, namely and Xinij, are still rather transparent. When the "spins" of the 
electron and the hole making the exciton are used - instead of the total kinetic momentum 
of the exciton -, these two parameters appear as a product of a spin part and an orbital 
part which have both a very simple physical meaning (see Eqs (22) and (25) for and 
Eqs (31) and (33) for \i ni j). The first one corresponds to all Coulomb interactions 
between the or (l,n) excitons when, on both sides, these excitons are made with the 
same electron-hole pairs (e,h) and (e',h'). The second one \i n ij has a purely fermionic 
origin and is simply related to the fact that the and (l,n) excitons can be made 
with different pairs, (e,h) (e',h') and (e,h') (e',h). This commutation technique allows 
to easily calculate any matrix elements between exciton states in an exact way. In this 
paper we have explicitly calculated the matrix elements of H between two-exciton states 
and we have shown that they cannot be related to the scattering rate of two excitons as 
previously proposed. 
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APPENDIX A : Calculation of the ^(q) 



Using Eq (17), 7a(q) reads 

7i*(q) = Putui(oiha x q) - Pu m (-a e a x q) , (68) 

with (5 vv i{w) defined by 

^(u) = (x v \e™- r,ax M = /£„(-u) , (69) 

and a x chosen to be the 3D Bohr radius. 

From Eq (69), we see that (3 vv i(0) = 5 uv i, the value u = being obtained either for 
q = or for a e> h = 0, i. e. rrih (or m e ) infinite. 

We can note that, if the v and v' states have the same parity, p uu >(u) = /?„„/(— u), so 
that 7/j(q) = for m e = m h : There is no direct scattering towards a same parity state 
if the electron and hole masses are equal. 

We can also note that the scattering 7^(q) depends on q = |q| only if the states 
are S states. 

Let us now calculate some values of these 7n(q), for 3D and 2D systems, when the 

states are S states. 
1) 3D case 

In 3D, the exciton relative motion wave function reads 

(r\x u ) = a?/ 2 <p<? D \p,e,<p), (70) 
where p = r/a x . The wave functions of the lowest energy S states are given by 

<fit D ] = e"7v^ , ^s D) = (2 - p)e^ /2 /v / 32T . (71) 

In 3D, Eq (69) reads for S states, 

&{u) = [ +0 °p 2 dp Pdhp rsm0d6e™» c ° se ^ D >(p,6,ip) ^, D \p,6,ip) . (72) 
Jo Jo Jo 

By inserting Eq (71) into Eq (72), we find 

022(«) = r, i > ( 73 ) 



(l+M 2 /4) 2 ' 

:i + m 2 



- 217/2 ^ (75) 

ft.,* W " 3 6 (1+4M 2 /9) 3 • ^ 
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From these /3^?\u), we can easily deduce the 7}f D ^(q) according to Eq (68), and the 
direct scattering according to Eq (23). 
2) 2D case 

In 2D, the exciton relative motion wave function reads 



(rW=«;Vi 2fl) (M). 



The lowest energy S state wave functions are given by 



(76) 



(2D) 
<Pls = 



2 3 / 2 /v^' 



-2p 



(2D) 



(2/3) 3/2 /v^l i 1 ~ 4 P/3) e" 2p/3 . (77) 



For 2D S states, Eq (69) reads 



2tt 



JO JO 



(78) 



Inserting Eq (77) into Eq (78), we find 

,(21?), - I 



)822(«) 



(1 + M 2 /16) 3 / 2 ' 

1 - (27m 2 /16) + (81m 4 /256) 



(79) 
(80) 
(81) 



(l + 9w 2 /16) 7 /2 

R[ 2D )( , 3 7 / 2 ^ 2 
PuM u ) 2 9 (l + 9« 2 /64) 5 / 2 ' 

From these /3^ u ?\u), we can deduce the Jif D \q) according to Eq (68) and the direct 
scattering according to Eq (23). For direct scatterings within the IS states only, we 
recover the result given in ref. (11). 
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APPENDIX B : Expressions of \i nij and in r space 

In order to show the equivalence between Eqs (24) and (25), it is simpler to start from 
Eq (25). Let us consider the V ee > term. By taking the Fourier transform of V ee >, namely 

Ke^E^V" 1 ^"^ > ( 82 ) 
q^O 

and the Fourier transform of the relative motion part of the wave function given in Eq 
(26), namely 




with = kj + a e Qi and K^ 1 = — kj + a^Qj, we can rewrite the V ee < term of Eq (25) as 
it.' (Ij) = 2 {xu, |kj) (zjk n ) (kija;^) (k^.) ]T Vq-^ ^ dr e d r/l dr e , dr h , 

k.i,k n ,k.i,kj q^O 

xe 4(-K-+K-+q).r e + (-Kf+K'').r h +(-K-+K j e -q).r e , + (-K^+K h ).r h ,] _ ^ 

The integral being equal to 5iq,Kf+q ^k^.k' 1 ^K«,K^-q ^kJj.k'm we immediately recover the 
first term of Eq (24). 

By transforming in the same way the other terms in Vhh', V e h>, Ve'h, it is easy to derive 
Eq (24) from Eq (25). 

In order to show the equivalence of Eqs (32) and (33), we can proceed similarly : By 
inserting Eq (83) into Eq (33) and by performing the integral, we immediately obtain Eq 
(32). 

We do the same to show the equivalence of Eqs (42) and (46) : We insert Eqs (82,83) 
into Eq (46) and perform the integral. This immediately gives the first term of Eq (42). 
The same transformation of the terms in Vhh', Veh', V e 'h leads to the three other terms of 
Eq (42). 
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APPENDIX C : Explicit calculation of F lsls i s i s (0, 0) 



From Eq (36), we get 



iW s i s (0,0)=]TKk|x ls >| 4 . (85) 

k 

The Fourier transforms of the relative motion wave functions (k|xi s ) are respectively 
given, in 3D and 2D, by 

V ' ' W(l + a? x k 2 ) 2 y ' 



2n a,, 



a x being the 3D exciton Bohr radius. 
In 3D we thus obtain 



(3D) 



rv rv\ 2 11 al x 2 dx 337Ta 3 

isisisuyO, 0) - — J q {1+x2)8 - ■ (88) 

In 2D we get 

FSSuu^ 0) = ^ r = ^# . (89) 
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FIGURE CAPTIONS 



Fig. (1) 

Direct Coulomb scattering of an i exciton into an I exciton, while a j exciton is 
scattered into an n exciton. As Coulomb interaction conserves the angular momenta of 
electrons and holes, we must have sj = Sj, s n = Sj, mi = to*, m n = m,j. As it also 
conserves momenta, we must have Qj + Qj = Q2 + Q„, between the centers of mass of the 
"in" and "out" excitons. Finally this direct Coulomb scattering contains the quantities 
7a(q) and j n j(— q), with q = Q2 — Q«, which characterize the scattering of a i/j exciton 
into a ^ state and a z/j exciton into a z/ n state, under a q excitation (see Eq (17)). 
Fig. (2) 

Bare exchange coefficient : The % and I excitons are made with the same electron but 
different holes ; we do have si = Sj but mi = rrij. 
Fig. (3) 

The dimensionless exchange coefficient A(p), given by Eq (38), as a function of P = 
p\2D, for 2D excitons and three different values of m e /mh- Thick solid line : m e /mh = 0. 
Thin solid line : m e /mh = 0.5. Dashed line : m e /mh = 1. 
Fig. (4) 

Right exchange Coulomb scattering (a), and left exchange Coulomb scattering (b). In 
the right exchange scattering, the exchange of holes making the excitons is made 
before the direct Coulomb interaction which then scatters the (p, r) excitons into (l,n) 
states. Note that the electron-hole interactions are between the (I, n) excitons but inside 
the excitons. Exchange and Coulomb processes conserving spins, we do have si = Sj, 
s n = Sj along with mi = rrij, m n = mi, the holes being crossed in the process. Note that 
the (I <-> n) change appearing in the definition of - and \i ni j generates a similar term 
in which the electrons are crossed instead of the holes. This restores the electron-hole 
symmetry. 
Fig. (5) 

The dimensionless difference between left and right exchange scatterings Ag(p), given 
by Eq (57), as a function of P = pA 2 D, for 2D excitons and three different values of 
m e /mh- Thick solid line : m e /mh = 0. Thin solid line : m e /mh = 0.5. Dashed line : 
m e /m h = I. 
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